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Abstract 

This paper presents an improved k-e model for low Reynolds number turbulence 
near a wall. The work is twofold: In the first part, the near-wall asymptotic behavior of 
the eddy viscosity and the pressure transport term in the turbulent kinetic energy equation 
axe analyzed. Based on these analyses, a modified eddy viscosity model with the correct 
near-wall behavior is suggested, and a model for the pressure transport term in the k- 
equation is proposed. In addition, a modeled dissipation rate equation is reformulated, 
and a boundary condition for the dissipation rate is suggested. In the second part of the 
work, one of the deficiencies of the existing k — e models, namely, the wall distance (e.g., 
y + ) dependency of the equations and the damping functions, is examined. An improved 
model that does not depend on any wall distance is introduced. Fully developed turbulent 
channel flows and turbulent boundary layers over a flat plate axe studied as validations for 
the proposed new models. Numerical results obtained from the present and other previous 
k-e models are compared with data from direct numerical simulation. The results show 
that the present k-e model, with added robustness, performs as well as or better than other 
existing models in predicting the behavior of near- wall turbulence. 

1. Introduction 

The k-e model is one of the most widely used turbulence models in engineering 
applications. Patel et al.M recently reviewed existing two-equation models that can be 
integrated directly to the wall. One of their conclusions was that the damping functions 
used in turbulence models, especially the one for the eddy viscosity, need to be further 
modified in order to improve model performance. In fact, as we shall see later, many 
existing k-e models do not provide the correct near- wall behavior of the eddy viscosity. 

Shih^ recently proposed a new near- wall k — e model based on asymptotic analysis. 
The present paper is a direct extension of that work. 

In the present paper, we will first analyze, in section 2, the near-wall asymptotic 
behavior of the eddy viscosity and the pressure transport term in the ^-equation, and 
in sections 3 and 4, propose models according to their near-wall behaviors. The model 
equation for the dissipation rate is reformulated following an argument similar to that of 
Lnmley,^ and a boundary condition for e is suggested. 

An asymptotic analysis shows that, in the near wall region, while the pressure 
transport term in the turbulent kinetic energy equation is small compared to the dissipation 
and molecular diffusion terms, it is much larger than the turbulent transport term, and 
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it plays an important role in the balance between the dissipation and molecular diffusion 
terms. This near-wall behavior is also observed in direct numerical simulation of fully 
developed channel flows (Mansour et alJ 32 ^ Kim et alJ 4 )). However, in existing k-e models, 
this pressure transport term is either ignored or lumped into the turbulent transport model. 
The present work introduces a model for the pressure transport term explicitly. 

Most of the existing k-e models for neax-wall turbulence use y+ (defined as u r yjv , 
where u T the friction velocity) as a parameter in constructing damping functions, with the 
Jones-Launder model being the only exception. While the use of y+ is perfectly fine for 
simple attached boundary layer flows, it is inconvenient in more general applications such 
as separated flows and flows with corners, where y + is not well defined. The Jones-Launder 
model has the advantage of avoiding y + ; however, the model is known to perform poorly 
in predicting near wall turbulent quantities, especially the turbulence kinetic energy. In 
the present work, a new damping function is derived based on asymptotic analysis (Section 
5). The new function is constructed upon a non-dimensional quantity that is independent 
of the coordinate system. 

The new models proposed in this paper were validated using direct numerical simu- 
lation data for fully developed turbulent channel flows and turbulent boundary layers over 
a flat plate. These nurperical results are reported in Section 6. Comparisons are also made 
with other popular k-e models implemented in the same computer code. The numerical 
results show that the present model, in general, performs better than the existing models 
while providing added robustness. 

2. Asymptotic Analysis 

To analyze the near-wall asymptotic behavior of the eddy viscosity and other tur- 
bulent quantities, we expand the fluctuating velocities and pressure in Taylor series about 
the wall distance as follows: 


u i = t>iy + C\y 2 -f djy 3 + . . . 
u 2 = c 2 y~ + d 2 y 3 + . . . 
u 3 = b 3 y + c 3 y 2 + <f 3 y 3 + . . . 
p = a p + b p y + c p y 2 + d p y 3 + . . . 


where the coefficients a p ,6i,c 2 ,... are functions of x,z and t . Using the continuity and 
momentum equations, Mansour et alJ 3 l derived the following relations between the coeffi- 
cients, 

2c 2 = —(^i,i + & 3f3 ) 

a pA = 2uc ! ( 2 ) 

a P)3 = 2 z/c 3 

The eddy viscosity is usually defined as 


(ie,Ujf) — + Uj,i) 



( 3 ) 
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where { ) stands for ensemble average and k = {v i u l )/2 is the turbulent kinetic energy. 
For plane shear flows, we can write from Eq. (3) 


_ ~{uv) 
dU/dy 


(4) 


and using Eq. (1), we obtain the near-wall asymptotic behavior of the eddy viscosity: 

u T = zihszly* + -M + c \c?) ±g^ifa)<£i) v 4 + 0(y*) (5) 

(hi) (&i ) 

That is, near the wall ut is 0(y 3 ). A correct eddy viscosity model should have this near- 
wall behavior. We shall see later that many existing models do not have this near-wall 
behavior. For later use, let us examine also the near-wall asymptotic behavior of the 
turbulent kinetic energy k and its dissipation rate e == viuijUij)- Using Eq. (1), we 
obtain the following relations for the k and e: 

k — - 1 - 2 " 3 -y‘ + + {^3^3))y 3 + 0(y 4 ) {6) 

~ = (b \ } + { 63 ) + 4((6}Ci) + {63 c 3 ))y + 0(y 2 ) (7) 

v 

In addition, the pressure transport term in the ^-equation, II = — becomes (using 

Eq.s (1) and (2)) 

n - — 2*/((6jC 1 ) + (6 3 c 3 ))y + 0(y 2 ) (S) 

The turbulent transport term in the fc-equation, — can be estimated as 0(y 3 ). 

Therefore, the pressure transport term is much larger than the turbulent transport term 
near the wall. 


3. Eddy viscosity model 

In this section, we will propose a model for the eddy viscosity using its near-wall 
behavior described in the previous section. In general, the eddy viscosity model can be 
written as 

vt = c u d (9) 

where v! and £! axe the turbulent characteristic velocity and length scale, respectively. 
Depending on how the velocity and length scales are specified, the eddy viscosity model 
can be a mixing length model, a one-equation ( k ) model or a two-equation (e.g. k-e) 
model. For example, in plane shear flows, Prandtl’s mixing length model specifies the 
characteristic velocity with ddUjdy. For near wall turbulence, the Van Driest mixing 
length model further damps the length scale to y [1 — exp(— y + /A)]. For more advanced 
mixing length models, see Baldwin and Lomax! 6 ], and King! 7 !. One-equation ( k ) models 
use k 1 ! 2 as the characteristic velocity, which is determined by the turbulent kinetic energy 
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equation. In two-equation models, e.g. k-e models, the length scale is usually specified 
by k 3 / 2 /e, where e is determined by a dissipation rate equation. In this paper we will 
concentrate on two-equation models, wherein the eddy viscosity is usually written as 


V T = 


1-2 


( 10 ) 


where = 0.09, and f M is a damping function. The form of the damping function 
is critical in such formulations, since the prediction of the mean velocity field depends 
primarily on the eddy viscosity model. Thus it is important for an eddy viscosity model 
to have the proper near- wall behavior. We have examined the near- wall behavior of eddy 
viscosity models based on various k-e model equations. The results are listed in Table 1, 
which shows that some of the k-e models do not have the correct near- wall behavior of the 
eddy viscosity, namely, = 0(y 3 ). 

The quantity k 3/2 /e is usually considered a characteristic length scale, (or the 
size) of the energy containing eddies. One expects that near the wall the size of these 
eddies to be of the order of the distance from the wall, O(y). However, Eq.s (6) and (7) 
show that k 3 / 2 /e is 0(y 3 \. Hence, k 3 / 2 /e is not an appropriate quantity to represent the 
length scale of the large eddies near the wall. We therefore introduce a new variable e: 


e 


e — v 


dk/dx, dk/dxi 

2k 


( 11 ) 


which has the following property: e approaches e away from the wall and is 0(y 2 ) near the 
wall. Therefore, k 3 / 2 /e is a proper quantity to characterize the length scale of the large 
eddies. With this length scale, the eddy viscosity should be written as 


V T = Cufu — 


( 12 ) 


Now in order for vy to have the correct near-wall behavior, the damping function f ^ 
must be 0(y) near the wall and approaches 1 away from the wall. The damping functions 
used in various k-e models are listed in Table 2. If we consider the presence of the wall 
as the main effect on the eddy viscosity, then we may assume / p is mainly a function of 
y + . The form of can be determined quite accurately if we know i/y, k and e from, for 
example, the direct numerical simulation. One may optimize the following simple form by 
numerical experiments: 


fn = 1 - exp(-aiy + - a 2 y +2 - a 3 y +3 - a 4 y+ 4 ) (13) 

The optimal values for channel flows are aj = 6 x 10 _3 ,O2 = 4 x 10~ 4 ,a3 = —2.5 x 
10 -6 ,a 4 = 4 x 10 -9 . It can be shown that this form of damping function does provide 
the required near-wall behavior. As will be shown later, the above constants are valid for 
general boundary layer flows. 


90 



4. Modeled k-e equation 

To complete the eddy viscosity model, we need the modeled equations for the tur- 
bulent kinetic energy and its dissipation rate. In this section we will analyze the near-wall 
behavior of the Jt-equation and propose a model for the pressure transport term with a 
proper near-wall behavior. The equation for the dissipation rate is also reformulated with 
a formal invariant analysis. 


4.1 ^-equation 

We start with the equation for the turbulent kinetic energy, 

k.t + Ujkj = D v + T + n + P-e (14) 

where D v , T and II represent the transport of the turbulent kinetic energy due to the 
viscosity, turbulent velocity and pressure, respectively. P and c are the rate of production 
and dissipation of the turbulent kinetic energy. The terms on the right hand side of Eq. 
(14) are defined as follows: 

D„ = vk tJJ 
T = —(knj)j 


11 = 


(15) 


P — (uiUj)Uij 
e = is(u,jUij) 

Using Eq.s (1) and (2), we obtain the budget of the /.--equation near the wall, 

— — = 0(y 3 ) 

Dt K ’ 

D„ = u({b]) + (bj)) + 6v((biCi) + {b 3 c 3 ))y -f 0(y 2 ) 

T = 0(y 3 ) (16) 

n = — 2i/((6ic 1 ) -1- {b 3 c 3 ))y + 0(y 2 ) 

P = 0(y 3 ) 

e = i/({6 2 ) + (6 2 )) +4z/((6 1 c 1 ) -f {b 3 c 3 ))y + 0(y 2 ) 

This budget shows that the term II is much larger than the term T, and that II cannot 
be neglected if we want the fc-equation be balanced in the near- wall region. However, the 
existing models either do not consider this term or simply combine it with the term T and 
model them as 


— (kuj)j — 


v JLy t . 

Ok J 


(17) 




In this paper, we propose a model for II which has a similar form to that of the standard 
turbulent transport model, but with a coefficient to ensure its correct near -wall behavior, 
Eq. (8). The proposed model form of II is 

Co 


n = 


VT 


/(i[l -exp(-y+)] a k 


*•'} 


(18) 
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where Co 0.05 is a model constant. In some existing k-e models, it is assumed that e = 0 
at the wall. In that case, in order to balance the term D v , a nonzero artificial term D 
must be added to the ^.'-equation. The form of D for various k-e models is listed in Table 
3. Finally, the modeled ^-equation becomes 


dk dk 

dt 3 dxj dxj 


9 ,/ , "Ty, M l TT ,dUi DU, dU x 

l(,, + 7- k ^ ]+n+ ^gr J + r I ■ ) a^~‘ +D 


(19) 


In the present model, D = 0, since e is nonzero at the wall. (The boundary condition for 
e will be discussed later.) 


4.2 6-equation. 


The exact dissipation rate equation is 


+ Vjtj = Dl + T< + Il< + PD (20) 

where D u> T and II* represent the diffusion rate of the dissipation rate due to the viscosity, 
turbulent velocity and pressure, respectively, and PD stands for the entire mechanism of 
the production and destruction of the dissipation rate e. The terms on the right hand 
sxle of the above equation are as follows: 


D i = u€ <jj 

T ( = -u(u, ik u lik uj)j 

n = ~{p,k u j,k),j (21) 

PD = — 2i/((u li r.u Jj t) -f {u kti u k j))Uij — 2v{v.jU lk )U lk] 

- 2v(u t}k u j ' k u iJ ) - 2 u 2 {i ii, k ju iikj ) 

The term II is usually neglected and the term T f is modeled as 


T* = 


l/'f 


-e.v 


( 22 ) 


To model PD, we define ^ by 


PD = -- t ** 
k 


At the level of the k-e model, we assume S' is a function of i/, u T , k, c, e, Uij and U t •*. 
Since S' is an invariant, it must be a function of the invariants that can be constructed 
from these quantities. Therefore we can write 


S' = S <(R t , 


vrfhjUij k 

. , wxU x% j k Ui t j k ^_) 


where R t is the turbulent Reynolds number k 2 /ue. We expand S' in a Taylor series about 
U,,j /( and vi/'pb ij k Uij k k j ee, and take only the linear terms. We obtain 


V"rU- U ■ l- 

* = >,J + rfavi 'rUijkUijk-p 


(23) 
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where the coefficients t/>o, V"’ 1 and v ,; 2 are in general functions of R t . Finally, the modeled 
dissipation rate equation becomes 


de de 
di +U] dx } 


-£-[{y + — )-| 1 ] + Ci/i jvrUtjUij - C 2 f 2 j + E 
OX] cr € OX] k k 


(24) 


where, C\ and C 2 are the model constants, and f\ and fi axe functions of Rt. The term 
E in the present model comes from the last term in Eq. (23): 


E = ui*rU itjk Ui tik (25) 

where we have taken rp2 == —T The form of E and C \ , C 2 , j\ and /2 for various k-e models 
are listed in Tables 3 and 4. 


4,3 Boundary Condition for e. 

Many of the earlier k — e models use e = 0 as the boundary condition for the 
dissipation rate. It is now generally agreed that this is not the physically correct boundary 
condition. However, controversy still exists in what boundary condition should be used for 
the dissipation rate. In some calculations, de/d y = 0 is used, which clearly has no physical 
background. Most models use the second derivative of the turbulent kinetic energy at the 
boundary as the boundary condition for e, as listed in Table 1. This condition comes 
directly from the ^-equation and is physically correct; however, it makes the problem very 
stiff and thus put very stringent restrictions on the choice of initial profiles for k and e. 
If the initial profile for k is not given correctly, the second derivative of k can become 
negative and cause the solution procedure to diverge. 

We propose the following boundary condition based on the asymptotic analysis. At 
y = 0, it is obvious from eqs. (6) and (7) that 


e = 2u 



(26) 


This expression is exact at the wall, and it does not add stiffness to the solution procedure. 

5. Deficiencies and Improvements of Existing Models 
5.1 Damping functions 

One of the deficiencies of the existing near wall models is that most of the wall 
damping coefficients are functions of a wall coordinate, such as the y+. This types of 
damping function works well only in the cases of attached boundary layer flows with simple 
geometries where y + is well defined. For practical engineering problems with corner flows 
or separated flows, some ad hoc treatment to the damping function must be made. The 
same is also true for the length scale in an algebraic model. The only exception to the above 
is the Jones-Launder model in which the damping function is a function of R t = k 2 fue. 

Although the Jones-Launder model has the advantage of independent of y + , it 
is known that this model does not predict correctly the near-all turbulence, especially, 
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the near-wall turbulent kinetic energy is underpredicted. One natural option could be to 
modify the J-L model such that it would predict the correct near-wall turbulent quantities. 
However, one basic difficulty met in such attempts is that the k -equation is very stiff. The 
fact that we are putting k 2 back into the equation by using the parameter R t aggravates 
the situation. 

To avoid the above difficulty and yet still achieve the same end, we introduce the 
following parameter: 

„ t / 4 

Ru = — , ( 27 ) 

where U is the total velocity. (Note: U is the total velocity in a coordinate frame fixed to 
the solid boundary.) From the results given in Section 3, since e approaches a finite value 
at the wall, it is obvious from eq. (1) that R u = 0(y 4 ) near the wall. Similar to eq. (13) 
of Section 3, we write the damping function for the eddy viscosity as: 

- 1 - exp(~aiK /4 ~ a 2 R]/ 2 ~ asRu) (28) 

with ai = 5 x 10~ 3 , a 2 = 7 x 10“ 5 , and a 3 = 8 x 10“ 7 . One can easily verify that with this 
damping function, the eddy viscosity has the correct near wall behavior, i.e., v? = 0(y 3 ). 

One point worth mentioning is the wide applicability of the above damping func- 
tion. Though developed with Shih model (Section 3, 4) in mind, it can be used with any 
existing k — t model that uses a non- zero boundary condition for the dissipation. This new 
parameter i? u , unlike R u by no means affect the stability of the solution procedure. 

5.2 Pressure transport term. 

In order to remove the coordinate dependency of the fc-equation, we replace the 
pressure transport term given in eq. (IS) by the following expression: 


where the model constant is readjusted to Co = 0.01. 

5.3 The formulation of e. 

In order to obtain the correct wall behavior for the eddy viscosity, we have in- 
troduced 6 in eq. (11), Section 3. Theoretical analysis shows that e is always positive. 
However, in numerical calculations, the value of e may become negative or even oscillatory 
due to round of errors. ( Depending on the accuracy of the numerical procedure, this may 
or may not be the case.) Here, an alternative definition of l is suggested: 


n = 


Cp VT j 

l/u <7* 


e = (1 - exp(-R] /2 ))e 


(30) 


This expression has the same near wall behavior as eq. (11) but is less likely to cause 
numerical instability. 
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With the modifications suggested above, the model constants need slight adjust- 
ment. The constants used in the present modified model are as shown in Table 4. 

6. Numerical Testing 

Flows with self-similar solutions are particularly useful for accurate model testing, 
because their solutions are independent of initial conditions, and one does not need to 
choose carefully the initial conditions for the k and e. In this paper, we use a fully 
developed channel flow and a flat plate boundary layer for model testing. These flows are 
the simplest wall boimded turbulent shear flows with self-similar solutions. However, the 
complex features of the turbulence, for example, the effect of the wall on shear turbulence, 
are present. In the case of the channel flow, the k-e equations form a one-dimensional 
problem, numerical calculations are easy and accurate. Recently, the measurements! 8 ' 
confirmed the accuracy of the direct numerical simulation data.M These data are used for 
model validations. 

In legends of the figures presented at the end of the paper, which will be discussed 
in detail in the following paragraphs, the word “present” refers to results obtained using 
the model suggested in Section 4 together with the new damping function and t given in 
Section 5, while the label “Shih” refers to results obtained using strictly formulations given 
in Section 4. 

6.1 Fully developed turbulent channel flow 

Let h be the half width of the channel, u T the friction velocity and Re r the Reynolds 
number defined a s u T h/v . Let £/, v x and y be the non-dimensional quantities, normal- 

ized by u r ,u^, u\/h,v and /i, respectively. The modeled equations for the channel flow 
become 


du 1-y 

— - = Re r 

dy l + i/ 


d f 1 r , __ \ , d k , / df*J , 2 1 

^W 1+(1+C, ;r j ^ )+ ‘' T( ^ ) jr r 


— e = 0 


d , 1 u T de , e <dU ^ 2 1 

+ 77 ^' +Cl k l ' T ^~dy) ~ 


„ r el ,d 2 U . 2 1 


Rel 


(31) 

(32) 

0 (33) 


where 


k 2 

V'T — f (J. R&T ~ 


e = e — 


(ik)2 

l dy l 

2k Rc t 


f n — equation(13) 

6e 

C = 

fni 1 -exp(-y+)] 


Re T k 2 )2] 


(34) 
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The boundary conditions are simple. 

At the wall, 



o 

II 

•'k 

II 

& 


and at the center of the channel, 

f _ (f) 2 

2k Re T 

(35) 

dk de 



dy dy 

(36) 

The main results from different k-e models for a channel flow with Re r 

= 180 


axe plotted in figures 1 4. All the calculations are compared with the direct numerical 

simulation data. Figure 1 shows mean velocity profile, Figure 2 shows the turbulent kinetic 
energy, Figure 3 shows the turbulent shear stress distribution, and Figure 4 shows the 
Dissipation rate. From these numerical results we reach the following conclusions: The 
model of Jones and Launder! 10 ! (JL) underpredicts the mean velocity as well as the peak 
value of the turbulent kinetic energy. Chien’s model! 11 ! performs better than the JL model, 
but it overpredicts the mean velocity near the center of the channel as well as the turbulent 
kinetic energy. In these two models, e = 0 at the wall is used as the boundary condition, 
so the dissipation rate near the wall cannot be correctly predicted. Lam and Breinhorsd 12 ! 
use a nonzero boundary condition for ( and have made some improvement for the mean 
velocity and turbulent kinetic energy compared with the results of the JL model. However, 
the shear is much overpredicted, and the dissipation rate near the wall is not correct. The 
model of Nagano and Hishida [13) presents a very good prediction for the mean velocity and 
shear stress, while the peak value of k is underpredicted. Their main modification to the 
JL model is a change in the damping function / M and the form of E. A zero dissipation rate 
at the wall is used. The numerical results from the Shih model and the present modified 
Shih model show improvements in the prediction of all quantities, including the dissipation 
rate. 

6.2 Boundary layer flows 

The boundary layer equations and the corresponding k- and e-equations are solved 
using a conventional semi-implicit finite difference scheme. In this scheme, the coefficients 
for the convection terms are lagged one step in the x-direction, and the source terms in 
the k— and e— equation are linearized in such a way that stability is ensured. 

In the present study, a grid of 100 points in the y-direction is used. The grid is 
stretched linearly with Axjj+i/Ayj = 1.05. The grids expands in the y-direction according 
the the boundary layer growth rate. 

The results of this calculation are presented in Figure 5 through Figure 10. 

Figure 5 and 6 show the comparison of the wall shear stress from various models to 
experimental data and some DNS data. As shown in Figure 5, at low Reynolds number 
(based on momentum thickness), J-L model overpredicts the shear stress while Chien 
model and NH model underpredict the shear stress. One common character of these three 
models is that they all used zero boundary condition for the dissipation, and thus unable 
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to predict the correct turbulence near wall behavior. At high Reynolds number (Figure 
6) the JL model still overpredicts the wall shear stress, while the others seem to do a fare 

job. 

Figure 7 and 8 are the results for Rc e = 1410, and Figure 9 and 10 are the results 
for Reg = 7700. From these results one can reach similar conclusions as we did from the 
channel flow case: The JL model generally underpredicts the peak value of the turbulent 
kinetic energy while overpredicts it in the inertia layer; the model also underpredicts the 
mean velocity profiles. The LB model and the NH model also underpredict the peak value 
of the turbulent kinetic energy near the wall. The three models mentioned above either 
have zero boundary condition for e or do not have the correct order of magnitude for eddy 
viscosity. The Figures show that, in general, the present model performs better than the 

existing models. 

Conclusions 

From the model testing, we conclude that the present k -e model has made consid- 
erable improvement over previous k-e models according to the comparison with the direct 
numerical simulation data. We find that the improvement is mainly due to the modified 
eddy viscosity model and the model of the pressure transport term in the fc-c quation. The 
proposed dissipation rate equation also shows a better near-wall behavior than the previ- 
ous ones. The correct boundary condition for e also seem to play an important role in the 
accurate prediction of the turbulence near wall behavior. 
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Table 1 Eddy viscosity and boundary condition for e in various k-e models 



Model 

t/'f 

BC: 


JL 

0(y 3 ) 

0 


Reynolds 

<?(y 5 ) 

"0 


LB 

0(y 4 ) 

dyl 


Chien 

0(y 3 ) 

0 


NH 

<V) 

0 


Shih 

0(y 3 ) 

*'(•#) 


Present 

0(y 3 ) 



Table 2 Damping functions used in various 

k-t models 

Model 
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Table 3 Model terms in various A>e models 
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Table 4 

Model constants in various k- 

e models 



Model 

c„ 


Ci 

C 2 

Ok 



JL 

.09 
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2.0 

1.0 
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Figure 2. Turbulent kinetic energy for channel flow, Re T = 180. 
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Figure 3. Turbulent shear stress for channel flow, Rc r 
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Figure 4. Dissipation rate for channel flow, Re r = 1 





Figure 5. Wall shear stress coefficient, flat plate boundary layer flows, low Reynolds 
number data. 
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Figure 6. Wall shear stress coefficient, flat plate boundary layer flows, high Reynolds 
number data. 
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igure 9. Turbulent kinetic energy for flat plate boundary layer, Re 


Mean velocity profile (Re, = 7700 ) 
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igure 10. Mean velocity profile for flat plate boundary layer y Rc$ 
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Low Reynolds Number Two-Equation 
Modeling of Turbulent Flows 


V. Michelassi* T.-H. Shih 
Center for Modeling of Turbulence and Transition 
NASA Lewis Research Center 

March 28, 1991 


Abstract 

A new k — e turbulence model that accounts for viscous and wall effects is 
presented. The proposed formulation does not contain the local wall distance 
thereby making very simple the application to complex geometries. The for- 
mulation is based on an existing k — c model that proved to fit very well with 
the results of direct numerical simulation. The new form is compared with nine 
‘different two-equation models and with direct numerical simulation for a fully 
developed channel flow at Re = 3300. The simple flow configuration allows a 
comparison free from numerical inaccuracies. The computed results prove that 
few of the considered forms exhibit a satisfactory agreement with the channel 
flow data. The new model shows an improvement with respect to the existing 
formulations. 
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MODELING OF NEAR- WALL TURBULENCE 


T.H. Shih* 

Institute for Computational Mechanics in Propulsion 
Lewis Research Center 
Cleveland, Ohio 44135 


N.N. Man sour 

National Aeronautics and Space Administration 
Ames Research Center 
Moffett Field, California 94035 


ABSTRACT 

Rpvnol!!! P 7 Se f ^ “ im P roved k( model a second order closure model for low- 

Reynolds number turbulence near a wall. For the k-e model, a modified form of the eddy 

viscosity having correct asymptotic near- wall behavior is suggested, and a model for the pres 
sure diffusion term m the turbulent kinetic energy equation is proposed. For the second order 

neTwalT beh ’ ^ “S * *?**** m ° dds ^ the Reynolds - stress equations to have proper 
near- wall behavior. A dissipation rate equation for the turbulent kinetic energy is also^efor- 

FuHv dev J H P T° S i fl 1$ SatiSfy leali2ability and wil1 not Produce unph^ical behavior 
di ' y t d ped , Chanael flows "e used ^ model testing. The calculations are compared with 
direct numerical simulations. It is shown that the present models, both the k-e model and the 

° SUre m ° dd ’ Perf0nD WeU in Promoting the behavior of the near wall turbulence 
Significant improvements over previous models are obtained. 
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ADVANCES IN MODELING 
THE PRESSURE CORRELATION TERMS 
IN THE SECOND MOMENT EQUATIONS 

Tsan-Hsing Shih and Aamir Shabbir 
Institute for Computational Mechanics in Propulsion 
and Center for Modeling of Turbulence and Transition 
Lewis Research Center 
Cleveland, Ohio 44135 

and 

John L. Lumley 
Cornell University 
Ithaca, New York 14853 


ABSTRACT 

In developing turbulence models, different authors have proposed various model con- 
straints in an attempt to make the model equations more general (or universal). The most 
recent of these axe the realizability principle (Lumley 1978, Schumann 1977), the linearity 
principle (Pope 1983), the rapid distortion theory (Reynolds 1987) and the material indif- 
ference principle (Speziale 1983). In this paper we will discuss several issues concerning 
these principles and will pay special attention to the realizability principle raised by Lum- 
ley (1978). Realizability (defined as the requirement of non-negative energy and Schwarz’ 
inequality between any fluctuating quantities) is the basic physical and mathematical prin- 
ciple that any modeled equation should obey. Hence, it is the most universal, important 
and also the minimal requirement for a model equation to prevent it from producing un- 
physical results. In this paper we will describe in detail the principle of realizability, derive 
the realizability conditions for various turbulence models, and propose the model forms 
for the pressure correlation terms in the second moment equations. Detailed comparisons 
of various turbulence models (Launder et al. 1975, Craft et al. 1989, Zeman and Lumley 
1976, Shih and Lumley 1985 and one proposed here) with experiments and direct numeri- 
cal simulations will be presented. As a special case of turbulence, we will also discuss the 
two-dimensional two-component turbulence modeling. 


107 


